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Development of a Physiologically-Based Pharmacokinetic 
Model for Sirolimus: Predicting Bioavailability Based on 
Intestinal CYP3A Content 

C Emoto 1 , T Fukuda 1 2 , S Cox 1 , U Christians 3 and AA Vinks 1 24 

Sirolimus is an inhibitor of mammalian target of rapamycin (mTOR) and is increasingly being used in transplantation and cancer 
therapies. Sirolimus has low oral bioavailability and exhibits large pharmacokinetic variability. The underlying mechanisms for 
this variability have not been explored to a large extent. Sirolimus metabolism was characterized by in vitro intrinsic clearance 
estimation. Pathway contribution ranked from CYP3A4 > CYP3A5 > CYP2C8. With the well stirred and Q gut models sirolimus 
bioavailability was predicted at 1 5%. Interindividual differences in bioavailability could be attributed to variable intestinal CYP3A 
expression. The physiologically-based pharmacokinetics (PBPK) model developed in Simcyp predicted a high distribution of 
sirolimus into adipose tissue and another elimination pathway in addition to CYP-mediated metabolism. PBPK model predictive 
performance was acceptable with C max and area under the curve (AUC) estimates within 20% of observed data in a dose 
escalation study. The model also showed potential to assess the impact of hepatic impairment and drug-drug interaction (DDI) 
on sirolimus pharmacokinetics. 
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Sirolimus (Rapamune®) was approved in 1 999 in the USA and 
in 2001 in Europe as an immunosuppressive drug intended 
to prevent organ rejection in renal transplantation. Sirolimus 
is an inhibitor of mammalian target of rapamycin (mTOR), 
a kinase enzyme which is an essential controller of cellular 
and organism growth that integrates nutrient and hormonal 
signals, and regulates multiple cellular processes. Recent 
studies have indicated that mTOR signaling is an impor- 
tant pathway in human diseases including cancer, diabetes, 
obesity, and neurological disorders. 1 Despite the increasing 
use of mTOR inhibitors in transplantation 2 and as a prom- 
ising agent in cancer therapies, 3 there is an unmet clinical 
need to better understand the dose-concentration-response 
and adverse events relationships. Sirolimus exhibits low and 
variable oral bioavailability with large interindividual phar- 
macokinetic variability. 4 To control some of the variability, 
therapeutic monitoring of sirolimus concentrations is recom- 
mended. CYP3A4, CYP3A5, and CYP2C8 have been sug- 
gested as the major enzymes involved in sirolimus hepatic 
metabolism. 5 However, the quantitative contribution of each 
pathway has not been characterized. CYP3A also represents 
the predominant drug-metabolizing system in the small intes- 
tine 6 and plays an important role in the intestinal first-pass 
metabolism. In addition, flavin-containing monooxygenase 
(FMO) contribution was explored as it may contribute to 
A/-oxidation of the sirolimus piperidine ring in a similar fash- 
ion as the FMO-mediated piperidine A/-oxidation described 
for phencyclidine. 7 After a single oral dose of 14 C-labeled 
sirolimus in healthy volunteers, 90% of the radioactivity was 
recovered from feces with little radioactivity recovered in 



urine (2.2%). 4 Sirolimus is reported to be a substrate of P-gly- 
coprotein (P-gp), but the contribution of P-gp to the binary 
excretion in humans is still unknown. 8 P-gp is localized in the 
canalicular domain of the hepatocyte and facilitates excre- 
tion of exogenous organic cations and xenobiotics into bile. 9 
Several metabolites were also detectable in feces and urine 
(Rapamune product information; Pfizer/Wyeth-Ayerst). Blood 
samples drawn as trough concentration from kidney trans- 
plant patients during maintenance dosing contained 56% of 
sirolimus-derived metabolites. 10 

After a single oral dose in healthy adults, a sirolimus con- 
centration decreases slowly with a half-life of 81. 5 h. 11 Siro- 
limus total body clearance is 278ml/h/kg with an apparent 
volume of distribution of 22.8 l/kg. 11 This large distribution 
volume implies extensive distribution to organs and tissues, 
and contributes to the long half-life. Sirolimus is partitioned 
extensively into blood cells with a blood to plasma ratio of 
35. 6. 12 This ratio is independent of sirolimus concentration 
in kidney transplant patients and exhibits a large variability. 13 
Sirolimus is also highly lipophilic with a log P of 4.3 (http:// 
www.drugbank.ca/drugs/DB00877) and the free fraction in 
plasma is 8%. 4 

Physiologically-based pharmacokinetics (PBPK) is a 
modeling technique that provides insight into the underlying 
mechanisms related to drug absorption, distribution, metabo- 
lism and excretion. Fundamental to PBPK is the separation 
of information on the system (i.e., human body physiology) 
from that of the drug (e.g., physicochemical characteristics 
determining permeability through membranes, partitioning 
to tissues, binding to plasma proteins or affinities toward 
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certain metabolizing enzymes and transporter proteins) and 
the study design (e.g. dose, route and frequency of admin- 
istration, concomitant drugs, and food intake). PBPK mod- 
eling is increasingly used in pharmaceutical research and 
drug development. In a recent report, it was noted that the 
office of clinical pharmacology at the US Food and Drug 
Administration received 25 submissions containing PBPK 
analysis between June 2008 and December 201 1 . 14 In addi- 
tion, several software packages have become commercially 
available for human PBPK modeling. The aim of the study 
was to first identify responsible CYP and FMO isoforms and 
quantify their contributions to overall sirolimus in vitro intrinsic 
clearance. Based on these data, we developed a sirolimus 
PBPK model to simulate the effects of tissue partitioning, 
CYP3A activity in the gut and liver, hepatic impairment, and 
drug-drug interactions (DDI) on sirolimus bioavailability and 
pharmacokinetics. 



RESULTS 

Characterization of drug-metabolizing enzyme 
contribution to the overall sirolimus metabolism and 
bioavailability 

The sirolimus metabolic stability assay with recombinant drug- 
metabolizing enzymes, CYPs and FMOs, showed metabolic 
contribution of CYP3A4, CYP3A5, and CYP2C8 (Figure 1). 
The decline in sirolimus concentration as compared to con- 
trols was the largest when incubated with CYP3A4 micro- 
somes, followed by CYP3A5 and CYP2C8. The other CYP 
enzymes did not show any decline. FMOs did not show any 
catalytic activity under basic condition (pH = 9.5) and neu- 
tral condition (pH = 7.4). The estimated in vitro intrinsic clear- 
ance (CLJ values of CYP3A4, CYP3A5, and CYP2C8 were 
9.33, 3.96, and 0.25 ul/min/pmol CYP at 0.3 umol/l sirolimus, 
respectively. The CL int of CYP3A4 was 2.4-fold higher than 
that of CYP3A5. 
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Figure 1 Metabolic lability of sirolimus in recombinant P450 and FMO enzymes. Sirolimus (3 umol/l) was incubated with recombinant drug- 
metabolizing enzymes. P450 concentration was set at 50 pmol/ml. Protein concentration of recombinant FMO enzyme was set at 0.5 and 
1 .0mg/ml in pH = 9.5 glycine buffer (a) and pH = 7.4 potassium phosphate buffer (b), respectively. Data represent the mean of duplicate 
determinations. 
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F g , respec- 
were calculated as 0.96 ±0.03 



Fractions of dose that escape hepatic (F h ) and intestinal 
(F ) first-pass metabolism were estimated based on the well- 
stirred model 15 and the Q gut model, 16 respectively, using CL int 
obtained as described and CYP3A content values reported 
by von Richter. 17 Hepatic and intestinal CYP3A content 
was responsible for the magnitude of F h and 
tively (Figure 2). F h and F g 
and 0.16±0.10, respectively (mean ± SD from 15 subjects). 
Accordingly, the bioavailability was predicted to be 15 ±9% 
(mean ± SD from 15 subjects). 

Development of a PBPK model to predict sirolimus 
pharmacokinetics 

An initial PBPK model was developed using free intrin- 
sic clearance values of each CYP isoform, free fraction in 
plasma, blood to plasma ratio, permeability in Caco-2 cells, 
physicochemical parameters (such as molecular weight, 
pKa, and log P), and default values of tissue-plasma partition 
coefficients (K) implemented in Simcyp based on population 
data for a healthy subject. This initial model was uninforma- 
tive and showed large differences in predicted vs. observed 
concentration-time profiles. 
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Figure 2 Sirolimus predicted hepatic and intestinal bioavailabilities. 
Hepatic (F h , a) and intestinal (F b) availabilities of sirolimus were 
predicted based on the well-stirred model and the Q ut model, 
respectively, using total CYP3A content data of von Ricfiter et al. 
(2004), in which genotypes of CYP3A5*1/*3 (open circle) and *3/*3 
(closed circle) were investigated. 17 



As part of the model optimization, the K value for each 
tissue was subsequently tested at 10, 100, and a maximum 
value of 1 ,000 as an initial setting for screening purposes to 
identify the tissue most contributing to sirolimus volume of 
distribution. Model performance, as assessed by comparing 
observed vs. predicted pharmacokinetic profiles, improved 
when /C p was set at the maximum value for adipose tissue. 
The /C p setting did not have a noticeable effect for any of 
the other tissues such as brain, gut, kidney, liver, heart, and 
lung. Subsequent sensitivity analyses exploring K values of 
adipose tissue over a range from 250 to 1 ,000 indicated the 
best fit with a value of 1 ,000 especially to reflect the distribu- 
tion phase (Figure 3a). As sirolimus shows a high blood to 
plasma ratio of 35. 6 12 , the K value of 1 ,000 for the adipose 
tissue is converted into tissue-blood partition coefficient of 
28, which was calculated as the K value of 1 ,000 divided by 
blood to plasma ratio of 35.6. 

From this analysis, it also became clear that another 
elimination route such as biliary excretion and/or non-CYP/ 
FMO metabolism as an additional intrinsic clearance path- 
way (CL int other ) in the PBPK model has to be assumed to 
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Figure 3 Sensitivity analysis of adipose distribution and 
intrinsic clearance from another elimination route on sirolimus 
pharmacokinetics. Systemic blood concentration of sirolimus 
following an oral administration of a single 15mg dose of sirolimus 
was simulated (solid line) by the initial model. The adipose/plasma 
partition coefficient (K) was changed from 250 to 1 ,000 by 250 (a). 
The intrinsic clearance of another elimination route (CL int other ) was 
applied to introduce with a value of from 0 to 300 by 1 00 n |Jl/min/1 0 6 
cells. Dots represent observed mean concentrations from 22 healthy 
subjects after an oral administration of sirolimus at fasting. 39 
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adequately predict the sirolimus pharmacokinetic elimination 
phase profile. Incorporation of CL int other values of 100-200 ul/ 
min/10 6 cells improved the model fit, as assessed by mini- 
mization of the sum of weighted residual squares (Table 1 
and Figure 3b). The predictive performance of the model 
with CL int other values of 100 (Model 1) and 200 ul/min/10 6 
cells (Model 2) were tested using observed clinical data in 
a dose escalation scenario of another cohort. 11 As judging 
by the sum of weighted residual squares value, Model 1 had 
better predictive performance than Model 2 and provided 
estimates closer to observed C and AUC values at doses 

max 

ranging from 0.3 to 8.0mg/m 2 (comparable with 0.6-1 6 mg, 
Figure 4). 

On the basis of the results of the sensitively analysis and 
iterative model fitting, Model 1 was defined as the base siro- 
limus model for healthy adults and was used in further DDI 
and effects of hepatic impairment evaluations. In the evalu- 
ation of the pharmacokinetic interaction between sirolimus 
and diltiazem in healthy adults, median (25-75 percentiles) 
of AUC values for 180 healthy adult subjects after a single 
oral dose of sirolimus alone and in combination with diltiazem 
were predicted to be 550 (357-778) and 771 (558-1135), 
respectively, while observed AUCs were 701 (561-893) and 
1030 (827-1610), respectively. Observed results were based 
on recalculation with numbers obtained by reading out the 
reported points in Figure 2 described by Bottiger et a/. 18 As 
for the AUC ratios, the median, 25-75 percentiles, and the 
geometric mean were 1.47, 1.35-1.58, and 1.48, respec- 
tively, while observed AUC ratios were 1.44, 1.19-2.46, and 
1 .58. Using the cirrhosis module (Child-Pugh grade A, B, and 
C) implemented in Simcyp, Model 1 also allowed exploration 
of the AUC change due to hepatic impairment and compari- 
son with reported sirolimus AUCs in patients with cirrhosis 
(Figure 5b). 

DISCUSSION 

This is the first comprehensive study on the metabolism of 
the mTOR inhibitor, sirolimus using 20 recombinant drug- 
metabolizing enzymes, and also the first PBPK modeling 
study of sirolimus. The present study indicates that the 

Table 1 Effect of varying model parameters on C max and AUC t 



K p of (Ml/min/10 6 C max AUC t 



Model 


adipose 


hepatocytes) 


(ng/ml) 


(ng x h/ml) 


SWRS 


Initial 


1,000 


0 


75.3 (112) a 


914(144) 


1.50 




250 


0 


79.2(117) 


1 ,287 (203) 


2.30 




500 


0 


76.6(114) 


1,148 (181) 


1.89 




750 


0 


75.7(112) 


1,016 (160) 


1.63 


Model 1 


1,000 


100 


66.6 (98.8) 


712 (112) 


0.26 


Model 2 


1,000 


200 


60.2 (89.3) 


580 (91 .3) 


0.35 




1,000 


300 


55.2 (82.0) 


488 (76.8) 


1.25 


Observed 13 






67.4±22.8 


635 ±225 





AUC, area under the curve; CL, clearance; K p , tissue-plasma partition 

coefficients; SWRS, sum of weighted residual squares. 

a Fold error (%). b Data represents the mean ± SD from PK data after fasting/ 



CYP3A subfamily plays an important role in the metabolism 
of sirolimus in terms of metabolic capacity and abundance. 
Sirolimus is predominately metabolized by CYP3As with 
CL int of CYP3A4 > CYP3A5. This result is consistent in part 
with a previous report that CYP3A5 exhibited a lower intrin- 
sic clearance than CYP3A4 on in vitro hydroxyl metabolite 
formation of sirolimus. 19 CYP2C8 also showed metabolic 
activity to sirolimus, although the contribution was consid- 
ered to be small based on the CL. t value and the relative 
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Figure 4 Comparison of C max (a) and AUC t (b) of sirolimus in healthy 
adults between predicted m and observed values. A virtual dose 
escalation study was carried out at various doses of 0.3, 1 , 3, 5, and 
8mg/m 2 (comparable with 0.6, 2, 6, 10, and 16mg, respectively) 
for a representative healthy adult with Simcyp version 11. In the 
study, the pharmacokinetics of sirolimus was simulated by the 
PBPK models (Model 1 and Model 2, see Table 1). Model 1 (closed 
circle) and Model 2 (open circle) included CL int other of 100 and 200 
ul/min/10 6 cells, respectively. The predicted PK'parameters as C max 
(a) and AUC t (b) were compared with observed data (mean ± SD) 
from a dose escalation study with the same dosing regimen in adult 
healthy volunteers. 11 The solid line represents the line of unity and 
the dashed lines represent 80 or 120% of the unity (an area within 
20% error). 
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Figure 5 Prediction of pharmacokinetic parameters in case 
studies, for a drug-drug interaction (a) and in patients with hepatic 
impairment (b) by the developed base model, (a) Sirolimus blood 
concentrations in a healthy adult after an oral administration of 
10mg of sirolimus alone (solid line) or with coadministration of 
1 20 mg of diltiazem (dashed line). A virtual drug-drug interaction 
study was carried out with the model for 1 80 virtual healthy adults 
as identical as possible to those in the clinical study in terms of 
age and male-female ratio of subjects using Simcyp version 1 1 
and each line shows the mean of individual data. In the study, 
the concentration-time profiles of sirolimus were simulated by the 
PBPK models (Model 1, see Table 1 and Figure 4). (b) Effect 
of hepatic impairment on AUC of sirolimus after a single oral 
administration of a 1 5 mg dose of sirolimus. This pharmacokinetic 
study was carried out for individual virtual healthy adults, 
and virtual patients with cirrhosis Child-Pugh grade A, B, and C 
(CP-A, CP-B, and CP-C, respectively) implemented in Simcyp 
version 11 (100 subjects for each population). Obtained AUC 
values from the virtual study with the PBPK model (Model 1 ) were 
compared with observed AUC values (mean ± SD) in a published 
study. 3040 

abundance of CYP2C8 in the liver and the intestine. In an in 
vitro/ in vivo extrapolation, the CYP abundance plays a criti- 
cal role when the magnitude of overall metabolic clearance 
is considered. It is well-accepted that CYP3As are predom- 
inantly expressed in the liver and the small intestine. 620 
The abundance of CYP3A4 is 4.6-fold higher than that of 
CYP2C8 in human liver microsomes according to meta- 
analysis data in Caucasians. 20 In addition, protein expres- 
sion of CYP2C8 in the small intestine was less than the 
detection level while the average abundance of CYP3A4 
and CYP3A5 was estimated to 62 and 23 nmol/gut in 31 
intestines in Caucasians, respectively. 1621 



The analysis of F h and F g indicates that the large vari- 
ability in bioavailability of sirolimus could be explained in 
large part by variable expression of intestinal CYP3A due 
to the relatively high contribution of intestinal metabolisms 
reflected by Fg value of 0.16. With respect to CYP3A con- 
tent, von Richter et ai. reported no correlation between 
hepatic and intestinal expression within subjects. 17 Using 
their CYP3A content data from the same patients, bio- 
availability of sirolimus was predicted to be 15% which 
well corresponds with the reported clinical value of 14%. 4 
In general, CYP3A4 substrates with a hepatic intrinsic 
clearance of over 100ml/min/kg tend to exhibit a low bio- 
availability due to substantial intestinal metabolism. 22 As 
the predicted hepatic intrinsic clearance of sirolimus was 
800 ml/min/kg, sirolimus clearly meets this characteristic. In 
addition, Lampen et ai. reported an eightfold interindividual 
variability of the intestinal metabolite formation rate of siro- 
limus. 23 Further investigations into the relationship between 
hepatic and intestinal CYP3A in the same patient would 
facilitate more accurate prediction. 

Sirolimus is a highly lipophilic compound with a high blood 
to plasma ratio and small free fraction in plasma. 4 Sensitivity 
analysis in this study suggested that the distribution of siroli- 
mus to adipose tissue could be the most critical parameter to 
predict the pharmacokinetics profile of sirolimus in humans. 
A large volume of distribution of sirolimus might be caused 
by the distribution into blood cell and adipose tissue. Siroli- 
mus has a long half-life of 81. 5 h after a single oral dose in 
healthy volunteers. 11 In general, the adipose tissue is a stor- 
age compartment for lipophilic compounds with slow release 
from this storage site resulting in long elimination half-lives. 
Poulin et ai. reported that lipophilicity and unbound fraction 
in plasma are mechanistic determinants of in vivo adipose 
tissue-plasma partition. 24 In kidney transplant patients, the 
distribution volume of sirolimus correlated with body weight 
and body surface area, and it is possible that this correlation 
may relate to the lipophilicity leading to the distribution into 
fat tissue. 13 

The established PBPK model in this study shows the pos- 
sibility of another elimination route. The contribution of this 
alternate route to the overall sirolimus oral clearance was 
calculated to be 16% and needs to be addressed in a future 
study. There are several possible explanations such as bili- 
ary excretion and/or non-P450/FMO metabolism. Since the 
majority of radioactivity, which contains not only sirolimus 
metabolites but also sirolimus itself, was found in feces after 
a single oral dose of 14 C-labeled sirolimus in healthy adults 
(Rapamune product information; Pfizer/Wyeth-Ayerst), 
it would be reasonable to assume direct biliary excretion 
in the PBPK model. Sirolimus is reported to be a substrate 
of P-gp which is one of efflux transporters expressed in the 
hepatocyte canalicular domain and facilitates excretion of 
xenobiotics into bile. 8 ' 925 Thus, the binary excretion of siroli- 
mus would be a possible factor to explain its low bioavilability. 
There are several efflux transporters expressed in the cana- 
licular domain of hepatocytes such as multidrug resistance 
associated proteins, bile salt export pump, and breast cancer 
resistance protein. 9 Except for breast cancer resistance pro- 
tein, the transport capacities of these transporters for siroli- 
mus and its metabolites have not been fully characterized. 25 
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Another possible explanation would be degradation of siro- 
limus through hydrolysis or formation of its ring opened 
isomer. 826 

A DDI study plays an important role in the risk assess- 
ment for drug safety and effectiveness. Bottiger et al. 
has investigated the DDI between diltiazem and siroli- 
mus in healthy adults since diltiazem, a calcium channel 
antagonist, is used to control hypertension in renal trans- 
plant recipients. 18 CYP3As contribute to the metabolism 
of both sirolimus and diltiazem. 527 In their study, median 
of increased AUC ratio (25-75 percentiles) by combined 
administration with sirolimus was 1.44 (1 .19-2. 46). 18 The 
median of predicted AUC ratio by the PBPK model was 1 .48 
and close to the observed value. This illustrates the poten- 
tial of the PBPK model to predict likelihood and magnitude 
of changes in sirolimus exposure when combined with a 
CYP3A4 inhibitor. 

It has been well established that environmental and dis- 
ease factors can have an important impact on the expres- 
sion and/or activity of human P450 enzymes. 28 CYP3A 
content in the liver was decreased by 75% in patients with 
severe chronic liver disease. 29 In a sirolimus pharmacoki- 
netics study, patients with Child-Pugh grade A, B, and C 
hepatic impairment showed -40, -48, and -60% oral clear- 
ances compared with healthy adults. 30 The PBPK model 
as described in this study also allowed exploration of the 
influence of liver disease severity on changes in sirolimus 
exposure, although AUC values were slightly overestimated 
(Figure 5b). In addition to P450 expression levels, overall 
drug exposure can be affected by several other factors. 
Zimmerman et al. (2008) reported that albumin, triglycer- 
ides, and hematocrit are decreased in patients with severe 
hepatic impairment compared with healthy subjects. 30 
These factors may explain the slightly overestimated AUC 
in patients. Further studies would need to consider these 
physiological changes due to hepatic impairment. Thus, the 
model would provide new insights into interpretation of dif- 
ferent pharmacokinetic phenomena which lead to further 
investigations about the physiological difference between 
healthy and disease populations. 

In summary, a mechanism-based PBPK model of siroli- 
mus for healthy volunteers was developed in this study. The 
predictive performance of the model was good as evaluated 
by predicted AUC values being within ±20% of the observed 
clinical data in dose escalation scenarios. The model also 
allowed exploration of the influences of DDI with diltiazem 
and hepatic impairment on changes in sirolimus exposure. 
The model required some assumptions such as distribution 
to adipose tissue and an additional route of elimination and 
warrants a further study. In addition, an oral solution was 
used as the sirolimus formulation in the simulations. With the 
commercially available tablet, differences in solubility and 
dissolution rate may be additional factors in the absorption 
process. Nevertheless, the current PBPK modeling study 
provides us with mechanism-based insights into important 
contributing factors of observed interindividual differences in 
sirolimus pharmacokinetics. This prototype PBPK model will 
also be helpful to guide the need for conducting in vivo DDI 
studies and population specific studies such as pediatric and 
disease patients. 



METHODS 

Materials. Sirolimus was purchased from LC Laboratories 
(Woburn, MA). Zotarolimus was purchased from Molcan 
Corporation (Toronto, Canada). Potassium phosphate buffer 
(500 mmol/l, pH 7.4) and Tris buffer (500 mmol/l, pH 7.5), 
NADPH regenerating system were obtained from BD Bio- 
science (Woburn, MA). Glycine buffer (1 mol/l, pH 9.5) was 
purchased from TOCRIS bioscience (Bristol, UK). Other 
reagents used in this study were commercially available and 
of analytical grade. 

Recombinant enzymes (insect cell, oxidoreductase + b 5 
insect control, CYP1A2, CYP2A6, CYP2B6, CYP2C8, 
CYP2C9*1 (Arg144), CYP2C18, CYP2C19, CYP2D6*1 
(Val374), CYP2E1, CYP2J2, CYP3A4, CYP3A5, CYP4A11, 
CYP4F12, CYP4F2, CYP4F3A, CYP4F3B, FM01, FM03, 
and FM05) were purchased from BD Bioscience (San Jose, 
CA, USA). CYP1A2, CYP2C18, CYP2D6, and CYP4A1 1 
were expressed with oxidoreductase, but others were 
expressed with both reductase and cytochrome b 5 . 

Enzyme assay. A substrate depletion assay was performed 
according to the method of Obach (1999) with slight modi- 
fications. 15 Sirolimus was incubated with recombinant CYP 
microsomes in an appropriate buffer in a shaking water bath 
at 37 °C. The protein concentration was adjusted at 0.4mg/ml 
by adding the control microsomes on the assay for determi- 
nation of CL int . Regarding FMOs, incubations were conducted 
under two different pH conditions, basic condition (pH = 9.5) 
and neutral condition (pH = 7.4). Basic condition (pH = 9.5) 
is optimum for the FMO metabolic activity, but the degree of 
sirolimus degradation at this pH is higher than that in neutral 
condition (pH = 7.4). 31 Protein concentration at pH = 7.4 was 
set at twice the concentration of that at pH = 9.5 to enhance 
the FMO activity. Sample collection was at serial time points, 
where the 0 min time point indicated the time when the NADPH 
regenerating system and substrate was added into the reac- 
tion mixture for the P450 and FMO assay, respectively. Sam- 
ples were extracted with methanol/0.2 mol/l ZnS0 4 (80/20, 
v/v) including zotarolimus as an internal standard. 32 After the 
removal of protein by centrifugation at 21 ,200gfor 10 min at 4 
°C, the supernatant was subjected to high-performance liquid 
chromatographic separation with tandem mass spectrometric 
detection (LC-MS/MS). Data of CL int represent the average 
of triplicate determinations (Supplementary Data S2 online). 

LC-MS/MS analysis. Sirolimus was quantified by an inte- 
grated on-line solid phase extraction (SPE)-high perfor- 
mance liquid chromatography-tandem mass spectrometry 
(API-3000, AB SCIEX, Foster City, CA) system according to 
the method by Koal et al. with slight modifications. 32 On-line 
enrichment was performed by the POROS R1/20 column 
(2.1 mm x 30 mm, 20 urn particle size, Life Technologies) as 
a SPE column. As an enrichment eluent, 50% (v/v) methanol 
was used and its flow rate was 200 ul/min at a rinsing step. 
At zero time, an aliquot of 5 ul supernatant of each sam- 
ple was injected to the LC-MS/MS system. The SPE elution 
and analyte transfer to the HPLC column (Luna 5 urn Phe- 
nyl Hexyl, 2 mm x 50 mm, Phenomenex) were performed by 
switching the six-port valve at 2.0 min. Column temperature 
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was set at 35 °C. The configuration offers a fast SPE elution 
in a back-flush mode and an analyte transfer to the HPLC 
column by high organic content of the eluent (MeOH/H 2 0, 
97/3, 10 mmol/l ammonium acetate, 0.1% acetic acid) at the 
flow rate of 200 ul/min. During this analysis step, the flow rate 
of SPE eluent was increased to 1000 ul/min at 2.3 min, and 
back to 200 ul/min at 4.4 min. The valve was switched back to 
the start position at 4.8 min. 

The MS detection was performed in the multiple reaction 
monitoring mode with the positive ion mode of electrospray 
ionization. The transitions of m/z, namely m/z of the proton- 
ated/deprotonated molecular ion and predominant product 
ion, are as follows: sirolimus, 931.6-864.3; zotarolimus, 
983.6-934.4 (as an ammonium adduct ion for both). 

The analytical method is linear over a concentration range 
of 10 nmol/l to 3 umol/l. The intra- and inter-day coefficients 
of validation were 0.30-1 1 % and 1 .6-7.5%, respectively. The 
lower limit of quantification was 10 nmol/l. 

Calculation. The peak area ratios of sirolimus to an internal 
standard, zotarolimus, were used for quantification. The area 
ratios were converted to percentage sirolimus remaining, 
using the t = Omin peak area ratio value as 100%. The slope 
of the linear regression from log percentage remaining vs. 
incubation time relationships (-k e| ) was used in the conver- 
sion to in vitro T 1/2 , values by in vitro T V2 = 0.693//c e| . Conver- 
sion to in vitro CL int (in units of ml/min/pmol CYP) was done 
using the following: 



CL 



0.693 



/'nw'fror i/2 (min)xCYP conc (pmolCYP/ml) 



0) 



Oral bioavailability (F ral ) is defined as F oral = F a x F x F h , 
where F is the net fraction of dose absorbed from the intes- 

a 

tinal tract, F g and F h are fractions of dose that escape intes- 
tinal and hepatic first-pass metabolism, respectively. In this 
study, F a was assumed to be 1 , F h and F g of sirolimus were 
predicted based on the well-stirred model and the Q gut model, 
respectively, using the in vitro data generated in this study 
and individual hepatic and intestinal CYP3A contents data as 
reported by von Richter ef a/. 15-17 
F h was estimated by the following equation: 



Q h + (f uB xCL uH . nt ) 



(2) 



where Q h , f uB and CI_ uHint are hepatic blood flow (20.7 ml/ 
min/kg), free fraction of sirolimus in blood, and unbound 
intrinsic clearance in the liver, respectively. 33 In order to 
calculate free fraction of sirolimus in blood, free fraction of 
sirolimus in plasma (f p = 0.08) and blood to plasma ratio 
(35.6) were used. 412 

The unbound intrinsic clearance of sirolimus in the liver 
(CL uHint ) was calculated based on the recombinant CYP 
data. The nonspecific binding value of sirolimus to micro- 
somal protein was predicted to be 0.45 using log P (4.3, 
http://www.drugbank.ca/drugs/DB00877) and the micro- 
somal concentration (0.4mg/ml) according to the in silico 



method by Hallifax and Houston. 34 This value was used to 
convert the apparent sirolimus clearance into the unbound 
clearance. 



CL uHJnt = S(' SEF rhCYPi xC ^rhCYPiJnt X AbUndanCe CYPi,H) 

/=1 

xMPPGLx Liver Weight 



(3) 



where ISEF rhCYPi , CL urhCYPiint , Abundance CYPiH and MPPGL are 
the inter system extrapolation factor in order to correct for differ- 
ences in activity per unit enzyme between recombinant human 
CYP (rhCYP) and human microsomes, unbound intrinsic clear- 
ance in rhCYP, abundance of each CYP in human liver micro- 
somes, and milligram protein per gram of liver, respectively. 
ISEF values of CYP3A4/5 and CYP2C8 were set at 0.24 35 and 
1.0 (a default value in Simcyp version 11, CERTARA, USA), 
respectively. MPPGL and liver weight were 40mg/g of liver and 
21 .4g of liver/kg body weight, respectively 3637 
F was estimated by the following equation: 



F = - 



Q 



^gut + ((j.gut XCL, 



(4) 



where Q gut and CL uGint are the gut blood flow and unbound 
intrinsic clearance in the gut, respectively. The fraction 
unbound in the enterocytes (f ) was assumed to be 1 
according to the method of Yang et a/. 16 Q gut is described as 
a hybrid of both permeability through the enterocyte mem- 
brane and villous blood flow as follows. 



Q 



. Qyilli X ^ L perm 
^villi ^^~perm 



(5) 



Q vjllj and CL perm are an actual villous blood flow (18 l/h) and 
permeability clearance of sirolimus through the enterocyte, 
respectively. The CL perm is defined as a product of effective 
intestinal permeability (P eff ) and a surface area of the small 
intestine. P eff and CL perm were calculated by the following 
equations: 



log P eff = 0.4926 x log P a 



app,Caco-2 



-0.1454 



CL =P eff x2;rx radium x length 



(6) 



(7) 



where P.. and P o ro . have units of 10 -4 and 10 -6 cm/s, 

eff app,Caco-2 ' 

respectively. P appCaco _ 2 of sirolimus has been reported 
to be 3.1 x10 -6 cm/s. 37 Intestinal radius of 1.75cm and 
the length of 6 m were used in this study according to Yang 
et al. w 

In addition, the unbound intrinsic clearance of sirolimus 
in gut (CL uGint ) is calculated based on the recombinant CYP 
data: 



CL uG , int = £(ISEF rhCYPi xCL urhCYPiint x Abundance c 



(8) 



where Abundance CYPiG is abundance of each CYP in human 
gut. 
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For the calculation of intrinsic clearance, CYP3A4 
and CYP3A5 expression data by von Richter et al. were 
used as the abundance data. 1738 Hepatic and intestinal 
CYP concentration was scaled to the microsomal pro- 
tein level and the organ level, where the following factors 
were used: liver, 157mg homogenate protein/g of wet tis- 
sue and MPPGL as mentioned above; intestine, protein 
content of enterocyte homogenate (1.2mg/g homoge- 
nate) and weight of enterocytes (720 g when a total vol- 
ume of enterocytes of 720 ml and a density of 1 g/ml are 
assumed). 1738 

PBPK model simulations. A PBPK model was developed 
using Simcyp software version 11 (Supplementary Data 
S1 and S3 online). The healthy volunteer database was 
selected and the representative and individual value of vir- 
tual healthy adults were used for the simulation. Especially, 
the simulations in the DDI study were conducted in virtual 
subjects with demographics as close as possible to those 
reported in the original clinical study in terms of age and 
male-female ratio of subjects. 18 The well-stirred model was 
selected to derive the whole blood hepatic clearance. The 
advanced dissolution, absorption and metabolism model 
was used for the absorption process, where oral solution 
was selected as a formulation of sirolimus. The following 
assumptions were made: there was no active drug trans- 
porter into or out of the hepatocytes and enterocytes 
involved, therefore any kinetic parameters related to a trans- 
porter were not used in the model. The sensitivity analysis 
was performed in one parameter basis to assess the influ- 
ence on a sirolimus blood concentration profile at a time. 
Using in vitro data of sirolimus available, the PK simula- 
tion was performed based on the full PBPK model and then 
the predicted concentrations-time curve was compared 
with observed data. Sirolimus concentration-time profiles 
observed in a clinical trial was obtained from graphs in the 
previous publication by scanning with GetData Graph Digi- 
tizer ver. 2.25. Model fitting was evaluated by visual inspec- 
tion and the minimum value of sum of weighted residual 
squares. Inhibitor profile of diltiazem and cirrhosis popula- 
tion profiles (Child-Pugh score of A, B and C) embedded in 
Simcyp version 1 1 were used for the prediction of DDI and 
the effect of hepatic impairments on sirolimus PK. 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

y Sirolimus is a low oral bioavailability drug with 
large pharmacokinetic variability. Sirolimus is 
extensively metabolized but the relative contri- 
bution of the pathways involved is not known. 

WHAT QUESTION THIS STUDY ADDRESSED? 

y Sirolimus in vitro intrinsic clearance was charac- 
terized using 20 recombinant enzyme systems 
and a PBPK model was developed to explore 
factors explaining PK variability. 



WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

/ CYP3A4, CYP3A5, and CYP2C8 are the 
main contributors to sirolimus metabolism with 
CYP3A4 activity 2.4-fold higher than CYP3A5, 
and 37-fold higher than CYP2C8. Variability in 
bioavailability can be attributed to variability 
in intestinal CYP3A content. The PBPK model 
could well describe sirolimus pharmacokinet- 
ics in healthy subjects and allowed explora- 
tion of P450 mediated drug-drug interactions 
and effects of hepatic impairment on systemic 
exposure. 

HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

y The PBPK model provides important 
mechanism-based insight in sirolimus dispo- 
sition and factors contributing to observed, 
between and within patient PK variability. This 
will facilitate individualized dose adjustments 
based on therapeutic drug monitoring data. 
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